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solved numerically under the covariant instantaneous approximation with the kernels containing 
scalar confinement and one-gluon-exchange terms. The masses for the heavy diquarks and the 
doubly heavy baryons are obtained and the non-leptonic decay widths for the doubly heavy baryons 
emitting a pseudo-scalar meson are calculated within the model. 
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I. Introduction 

The past few years have seen many important developments concerning hadron colliders, especially the advent 
of the LHC. Recently, the discovery of and searches for double charm baryons have been reported by various 
experimental collaborations [ij-Ql ■ One is convinced that more and more doubly heavy baryons will be observed 
in the near future. Consequently, it is an urgent task for theorists to investigate the properties of these states. 

On the other hand, the existence of three valence quarks in a baryon makes the theoretical study much more 
complicated than the case of mesons. People have suggested the presence of diquark structure in a baryon and 
studied the properties of heavy baryons in such a picture . No matter whether the diquark is a real physical 
object or simply a theoretical approximation, this picture reduces the three body system to a two body problem 
which is much simpler for investigation. 

In recent years, heavy quark effective theory (HQET) has been widely used in the study of doubly heavy 



baryons 10l-ll7l|. Two heavy quarks are reasonably bound into a color-antitriplet heavy diquark whose radius 
is much smaller than the typical scale {1/Aqcd) of the nonperturbative QCD interactions in the heavy quark 
limit (toq S> Aqcd, tt^q denotes the heavy quark mass). The leftover light quark involved in the baryon moves 
in the color field induced by the heavy diquark. Unlike the heavy quark and light quark system, the internal 
motion in the heavy diquark can not be ignored even at leading order in the 1 / mq expansion. This is because 
the relative momentum in the heavy diquark is not simply 0{Kqcd), but ^ cij.'rnQ as calculated in the Coulomb 
potential model [3| • 

As a formally exact equation to describe the relativistic bound system, the Bethe-Salpeter (BS) equation 
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2Cl. However, it 



was initially formulated in Minkowski space based on the relativistic quantum theory 
is difficult to solve the BS equation in Minkowski space due to its singular behavior. In order to overcome 



* E-mail address: mhweng@mail.bnu.edu.cn 

t E-mail address: xhguo@bnu.edu.cn 

t E-mail address: anthony.thomas@adelaide.edu.au 



2 



this difficulty, with the so-called "Wick rotation", the formalism for the BS equation in Euclidean space was 
developed and investigated in detail 2ll424|. With the perturbation theory _integral representation, the BS 
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Recently, based on the 



equation was solved in Minkowski space for the scalar and fermion systems 
Nakanishi integral representation of the BS amplitude and the projection of the BS equation on the light- 
front plane, a new method for solving the BS equation in Minkowski space was proposed and was applied to 
study the electromagnetic form factor 31-3^. In another formalism to solve the BS equation, the covariant 



instantaneous approximation is adopted in the kernel. In recent decades, this formalism has been successfully 
used to investigate heavy mesons, heavy baryons, and exotic states 0, 

Based on the diquark picture of the composition of the doubly heavy baryon, we will establish the BS equations 
for both the heavy diquarks and the doubly heavy baryons in the leading order of a l/mg expansion. Motivated 
by the potential model, the kernel for the BS equation is assumed to be composed of the scalar confinement and 
one-gluon-exchange terms 45|. We will solve the BS equations numerically under the covariant instantaneous 
approximation 4l|-|4j|. Since the heavy diquark is not really a point object, a few form factors for the effective 



vertex of the heavy diquark coupling to the gluon are introduced to reflect the inner structure of the heavy 
diquark. These form factors will be expressed in terms of the BS wave functions obtained for the heavy diquarks. 
Finally, we will calculate the non-leptonic decay widths for the doubly heavy baryons emitting a pseudo-scalar 
meson in the BS formalism. 

The remainder of this paper is organized as follows. In Sec. II, we establish the BS equation for the heavy 
diquarks in the leading order of l/mg expansion. We also give the normalization conditions for the BS wave 
functions for the heavy diquarks in this section. In Sec. Ill, the form factors for the effective vertex of the 
heavy diquark coupling to the gluon are derived from the BS wave functions obtained for the heavy diquarks. 
In Sec. IV, we establish the BS equation for the doubly heavy baryons at leading order in the l/mg expansion. 
The normalization conditions for the BS wave functions for the doubly heavy baryons are also given in this 
section. In Sec. V, the non-leptonic decay widths for the doubly heavy baryons emitting a pseudo-scalar meson 
are calculated in the BS formalism. Sec. VI is reserved for our summary and some discussions. 



II. BS equation for heavy diquarks 

In general, the parity of a ground state baryon is positive. Since the parity of quark is supposed to be positive, 
the parity of the diquark involved in a ground state baryon should be positive. Due to the Pauli principle, two 
quarks with the same flavor can only constitute an axial-vector diquark. On the other hand, two quarks with 
different flavors can constitute either a scalar diquark or an axial- vector diquark. It can be easily shown that a 
heavy diquark which is in the ground state can not be a tensor diquark. 

Suppose two heavy quarks Qi and Q2 (with masses mg-^ and mg^ respectively) compose a ground state heavy 
diquark. Define two ratios Ai = mq^/ {itiq-^ + raq^) and A2 = "^Q2/("^(3i + ''tiQ2)- The BS wave function for 
the heavy diquark is defined as follows: 

XPo{^i.X2)c.p = e*^''=(0|r^i(a;i)j,V2(a;2)'^|-Pz3,fc) 

= e-^-^l ^xp.b)a,e-''^ (1) 

XPn{x2,Xi)pc. = £*^'^FD,fc|TV'2*(x2)^^Vi*(a:ira|0) 

= e^^-^l ^xp.b)/5ae^''^ (2) 
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where ■i/'i and V'2 stand for the field operators of the heavy quarks Qi and Q2, respectively, i, j, k represent 
the color indices, a and (3 represent the spin indices, X = AiXi + X2X2 is the coordinate of the heavy diquark 
mass center, x = xi — X2 is the relative coordinate of the two heavy quarks, Pd is the momentum of the heavy 
diquark, and p is the relative momentum between the two heavy quarks. 

The BS equation for the heavy diquark can be written in the following form (details can be found in Ref. 

): 

XPoip) = S{p,)(g>S{p2) J ^[r(E>j,K^''^ip-p')+I(E>IK^'f\p-p')]xPoip'), (3) 

where pi = XiPd + P and p2 ~ X2PD ~ P are the momenta of heavy quarks Qi and Q2, respectively, S{pi) 
and S{p2) are the propagators of heavy quarks Qi and Q2, respectively. K^^^^ and K'^'^^^ are the one-gluon- 
exchange and scalar confinement terms of the kernel for the BS equation given by (after imposing the covariant 



instantaneous approximation 4l|-|44|): 



and 



where as and k are the coupling parameters related to one-gluon-exchange and scalar confinement terms, 
respectively, pt is the transverse projection of the relative momentum (p) along the heavy diquark momentum 
(Pd) (see the definition below Eq. the second term of K^'^^^ is introduced to remove the infrared singularity 
near the point p'^ = pt, and the small parameter /x is introduced to avoid the divergence in the numerical 
calculations. This kernel is motivated by the potential model which has been successfully applied in mesons 



45 [ . Furthermore, we assume that the kernel of the heavy diquark is related to the meson by the one-half rule 

30- 



Eq. ([3]) can be written in a more usual matrix form as 

xlM = S{p2) J i]^[-7^x?,(y)7pi^(^^)(pt-p;)+xLb')^'^^Hft-p;)]^(-Pi), (6) 

where xpd (p) — ^XPd (p) the charge conjugation matrix) and the superscript T represents the transpose 
of the spinor indices. 

In the leading order of a 1/toq expansion, the heavy quark propagators {S{pi) and S{p2)) can be written as 

SiPi) = + ^, (7) 

2wQi [XimD +pi ~ ojQ^ + le) 

and 

Si.P.) = -^.^^^^^^^^^ -y (8) 

2^Q2 \-X2mr1 +Pi+ - le) 

where ttid and vd are the mass and velocity of the heavy diquark, respectively, pi = p-v^ andp^ = p'^ ~Piv'^ are 
the longitudinal and transverse projections of the relative momentum (p) along the heavy diquark momentum 



{Pd), respectively, the energy (jJq-^,^) = \/™Q "Pt, and e is the infinitesimal. As studied in Refs. id. 



1(2) - y '"51(2) 

\pt\ is 0{a^mQ). Consequently, \pt\/mQ terms should not be neglected in calculations carried out to leading 
order in the l/mq expansion. 
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Substituting Eqs. ([7]) (|H]) into Eq. ([SJ, one finds the following two constraints for the BS wave function for 
the heavy diquarks: 

^DXIJP) = XIJP)^ (9) 

xIAp)^d - -xIAp)- (10) 

Then, taking these constraints into account in the BS equation for the positive parity and zero angular 
momentum ground state of the heavy diquark system, the BS wave functions for the scalar and the axial-vector 
heavy diquarks can be parametrized in the following forms, respectively: 

xIAp) = i^D + lh'fu (11) 

and 

xPJip) = ii'D + l)t^f2, (12) 

(r) 

where is the r-th polarization vector of the axial-vector heavy diquark, /i and /2 are the Lorentz-scalar 
functions of , pi, and = mjj . 

After some algebra, we find that the BS scalar wave functions for both the scalar heavy diquark (/i) and the 
axial- vector heavy diquark (/2) satisfy the same integral equation as follows: 

HPt)- . "^^"^^ , r / M^[V^''Hpt~p't) + V^'fHpt--p'tmp't), (13) 



^^QiWQ2(-mD +a;Qi +wqJ y (27r)3 

where we define f{pt) = J ^/i(2)(p), V^^s) ^ and V^'=^^ = ~iK'^''f\ 

In general, the normalization condition for the heavy diquark can be written as (after imposing the covariant 
instantaneous approximation on the kernel) 



where Zi(2) and ji(2) represent the color indices of the heavy quarks, (5*^*^ = ^h^jl ~ ^h^h' ^'^"^ I^p^''^''^ {PiP') 
stands for the inverse of the four point function. 

Now, it is straightforward to obtain the normalization condition for the BS wave function for the heavy 
diquark as the following: 

~\j (^^{'^''[^^i'^(^^i)^^^o(^*)'^'^^^2)Xp^(pt)'5'(-Pi)^] 

+1r[\^S{-p,)^f^{pt)S{p2)iS{p2npM]^ 1' (16) 
where e = (1, 0), Xp^ (pt) and Xp^ (Pt) a-re the transverse projections of the BS wave function given by 

ipjpt) - -^siP2r'x^>Ap)s{^Plr\ (17) 

and 

XpliPt) = cf_pJ-pt)C-\ (18) 
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respectively. 

For the scalar heavy diquark, the transverse projections of the BS wave function {xpoiPt) XpjjiPt)) can 
be obtained from Eqs. ([TTl) (fT8|) as follows respectively: 

ipr, (Pt) = L iPth' + Is, {pt)^Dl\ (19) 



and 



where 



and 



xpjpt) = fsApth' ~ Lipt)^Di', (20) 



fsM) = I ^[v('^\pt-p't) + 'iv'^'^Hpt~PtmPt), (21) 



(22) 



After substituting Eqs. ([TO]) (PO)) into Eq. ([TC]). carrying out the trace calculation, and integrating out the 
longitudinal momentum pi , we obtain the normalization condition for the BS wave function for the scalar heavy 
diquark as the following: 



-[/s,(ft) + /s.(ft)]' = l, (23) 



where Ed ~ Pd ■ £■ 

For the axial-vector heavy diquark, the transverse projections of the BS wave function (xp^ {Pt) and Xp^ ipt)) 
can be obtained from Eqs. (IT71) as follows respectively: 

XpM) = f.APt)t^ + L{pt)fDt\ (24) 



and 



where 



xZ'\pt) = -LiPt)t^ + LAPt)fDe\ (25) 



and 



~fvAPt) = / ^[V-'^Kpt-p't) + 2V^^f\pt-p[)]f{p[), (26) 

L{Pt) ^ J ^V('^Hpt^Pt)fiPt)- (27) 
Analogously, the normalization condition for the BS wave function for the axial-vector heavy diquark is given 

by 

/ TTTTT — 5 2 2~( \ ; ^2 — [hAPtj + JvAPt)\ - t. 

In the numerical calculations, we take the constituent masses of the heavy quarks to be mi, = 4.88 GeV and 
rric = 1.486 GeV which were obtained by fitting the real spectra of charmonium and bottomonium in Ref. [4^. 
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The parameters in the kernel as = 0.4 and k = 0.18 are determined by fitting the experimental data for heavy 
meson spectra . In order to solve the integral equation ([T^ , we discretize the integration region into n pieces 
(with n sufficiently large). Then Eq. (|13p becomes an eigenvalue equation for the n dimensional vector /. After 
solving the eigenvalue equation, the heavy diquark masses are obtained and are displayed in Table ID We find 
that the heavy diquark masses are independent of the heavy diquark spin and only determined by the flavors 
of the constituent heavy quarks. In Fig. 1, the normalized BS scalar wave functions for the heavy diquarks are 
shown. It can be seen that the amplitudes of the BS scalar wave functions do not distinguish the different spins 
of the heavy diquarks. As discussed in Refs. 



IC 



unlike the heavy quark and light quark system, only the 
spin symmetry survives (when the diquark mass is below ~ 10 GeV) in the leading order l/mg expansion for 
the heavy diquark system. Our results are consistent with this statement. 



TABLE I: Values of the heavy diquark masses used here. 



niQ, {GeV) 


4.88 


4.88 


1.486 


mq^ {GeV) 


4.88 


1.486 


1.486 


mo {GeV) 


9.80 


6.55 


3.23 




0.5 1.0 1.5 

IPJ (GeV) 



2.0 



FIG. 1: The normalized BS scalar wave functions {f{pt)) for the heavy diquarks. The solid, dashed, and dotted lines 
are for the heavy diquarks composed of double b quarks, b and c quarks, and double c quarks, respectively. 



III. Form factor of heavy diquark coupling to gluon 



Since the heavy diquark is not really a point object and its radius is enhanced by In niQ with respect to 
1/toq, we introduce a few form factors for the effective vertex of the heavy diquark coupling to gluon to reflect 
the inner structure of the heavy diquark. ^ 

The effective current for the scalar heavy diquark coupling to a gluon is given as follows [5]: 

A'' 



(29) 
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o 
o 
o 



FIG. 2: The schematic diagram for the heavy diquark coupling to gluon. Qi(pi'') and Q2{P2^) stand for the heavy 
quarks Qi and Q2 with momenta p^'' and '1 respectively. 



where gs is the couphng constant of the strong interaction, A° (a = 1, 2, • • •, 8) denote the Gell-Mann color 
matrices, Fs{Q'^) is the form factor for the effective vertex, P^. and P'^^ are the momenta of the initial and 
final heavy diquarks, respectively, and is the square of the momentum transfer. 

On the other hand, the effective current for the scalar heavy diquark coupling to a gluon can be written as 
the following in the BS equation formalism (The Feynmann diagram for the vertex is shown in Fig. 2): 

= ^g,^{M^ + M^), (30) 



where 



xTr[5(-pi)^^i^ iPt)SiP2)^P„^ ip't)Si-p[)v% (31) 



and 



= --{2^)H\Pd- - PD,-k) ( f ^-KjC^tt 

2 Qy ' y 'J (27r)4 J (27r)4^ 



)'SHpi~p'() 



xTr[5(-pi)^^i^ iPt)S{p2)v''S{p'^)ip,^ {p'^)h (32) 

where k denotes the momentum carried by the gluon, p and p'^"-* in this section denote the relative momenta 
of final and initial heavy diquarks, respectively. 

Then, comparing Eq. ((29)) with Eq. (pO|. one can get the form factor (Fs(Q^)) in the leading order of a l/mg 
expansion as follows: 

(27r)4 f d^pt r rnlmg^iUipt) + fs2{pt)][fsi{p't) + ~fs2{p't)] 



3mn I {2nY L 



3m_D J (27r)3 I ujQ^ujQ^ujQ^{mD - ujqi - 0JQ2){mD - ojq^Q - - lcos9 - WgJ 

mQ,ml^[fsi{pt) + fs2{pt)][fsM) + fs2{p'l)] 



^Qi^Q2^Q2^'m'D - - ujQ^){'mD - ujq^lo + \ptW i^'^ - lcos6 

where the velocity transfer Cj — ■ vOf (fD; and are the velocities of the initial and final heavy diquarks 
respectively), pf"'' — p'^"^ ■ and p^"^ = p'^"^ —p^'^VBi are the longitudinal and transverse projections of 
the initial heavy diquark relative momenta {p''^"^) along their momenta [Poi) respectively, uj'q = \/"^q ~ P't^ 



Jq^ = Y^^Qj — p'P , and is the angle between pt and Pj^.^i— Pp,. — v^^Poi ■ vdj)- 
Now, let us turn to the effective current of the axial- vector heavy diquark coupling to gluon 

A" 



+ P'D.)Fv.m - (P'dX'^ + PB.9'')F.AQ') + PB.P^.iP^, + P'D.)Fv.m , (34) 
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where Fy^{Q'^), Fy^{Q'^), and Fy^{Q'^) are the form factors for the effective vertex. 

Analogously, the effective current of the axial- vector heavy diquark coupling to the gluon can be written as 
follows in the BS equation formalism: 

J"f/3 = ^g, — {M^^'^ + M^^'^), (35) 



where 
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xTr[5(-pi)^^i';(pO^(p2)^P,^(p;)^(-p'i)«^], (36) 



and 



xTr[5(-pi)^S"(P*)^b2)^'^5(pg)^^^^(p;')]- (37) 

As discussed in Ref. 5], the contribution of the Fv^{Q'^) term is suppressed at small and intermediate mo- 
mentum transfer, Q^, since such a term is multiplied by high powers of momenta. Consequently, F^^(Cp) is 
ignored in our calculation. Comparing Eq. ([M)) with Eq. ([35l) . we derive the other two form factors {Fy^iQ"^) 
and Fy^{Q'^)) in the leading order of a '^/mq expansion as follows: 

(2^)4 r d^pt r "lQi"^Q2 [fv^ [Pt) + fv, {Pt)] [fv, (P't) + fv2 (Pt)] 



12TOn ./ 27rP L 



12toi3 J (27r)3 ^ujQ^ojQ^ujQ^lmD - ujq^ - ujQ^){mD - loq^^ - \ptW Cj'^ - IcosO - WgJ 
'7^Ql"^L [/^i {Pt) + {pt)] [fv, (Pt) + fv2 (Pt)] 



Fy,{Q^) = 0. (38) 

It can be seen that the form factors for both the effective vertex of the scalar heavy diquark coupling to the 
gluon and the effective vertex of the axial- vector heavy diquark coupling to the gluon are equal to each other 
in the leading order of a l/mg expansion. So, we redefine the form factors as F{Q^) = Fs{Q^) — Fy^lQ"^) for 
convenience. 

It is well known that when — > 0, the heavy diquark is seen by the gluon as a point particle without 
the inner structure, and hence the form factor for the effective vertex should be normalized to unity. When 
— > cxD, the gluon can see the individual quarks inside the diquark, and hence the form factor for the effective 
vertex should approach to zero. We calculate the form factors with the BS wave functions obtained numerically 
for the heavy diquarks. The dependence of the form factors on the square of the momentum transfer is shown 
in Fig. 3. We can see that the behavior of our results coincides with the tendency of the above physical picture. 



IV. BS equation for doubly heavy baryons 



As discussed in Sec. I, the doubly heavy baryon can be regarded as a bound state of a heavy diquark and a 
light quark in the heavy quark limit. 

Let us define the ratios rji = nii/{mi + m^)) and •q2 — Tno/imi -I- mjj) {mi is the light quark mass) for 
convenience. The BS wave function for the doubly heavy baryon composed of a scalar heavy diquark and a 
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Q' (GeV') 



FIG. 3: The normalized form factors for the vertices of the heavy diquark couphng to gluon as a function of the square of 
the momentum transfer. The sohd, dashed, and dotted hues represent the heavy diquark composed of double b quarks, 
b and c quarks, and double c quarks, respectively. 



light quark is defined as the following: 

xp{yi,y2) = (o|V'i(yi)0D(y2)|f') 

= ^-^""/(§4Xp('Z)e--, (39) 

where ipi and stand for the light quark field and the scalar heavy diquark field, respectively, Y = rjiyi +7?22/2 
is the coordinate of the doubly heavy baryon mass center, y = yi — y2 is the relative coordinate, P is the 
momentum of the doubly heavy baryon, and q is the relative momentum between the heavy diquark and the 
light quark. 

It is straightforward to derive the following BS equation for the doubly heavy baryon containing a scalar 
heavy diquark and a light quark: 

Xpiq) = Siiq,) J ^GiP,q,q')xp{q')SDiq2), (40) 

where qi — r]iP + q and 52 = V^P ^ 1 a-re the momenta of the light quark and the heavy diquark, respectively, 
Si{qi) and SDiq2) are the propagators of the light quark and the heavy diquark, respectively, G is the kernel 
which is, motivated by the potential model, given by 45 1 

~ lG{P,q,q') = I ® IViiq^q') +-f^,®T^V2{q,q'), (41) 

where ~ {q!^ + q2^)F{Q'^) is the effective vertex of a gluon with two scalar heavy diquarks, which has been 
derived in Sec. Ill, Vi and V2 are the scalar confinement and one-gluon-exchange terms given in the following 



respectively (after imposing covariant instantaneous approximation 4l|-|44j): 



and 



[-{'It - q'tV + f^^V J {277^ [-{qt - hy + fi^Y 

V2{qt ~ q't) = -^7 ^5 J, 43 

3 [qt - q'tV - 
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where q^'^s the transverse projection of the relative momentum (q) along the baryon momentum (P), the second 
term of Vi is introduced to remove the infrared singularity near the point qt = q't, and the small parameter /i is 
introduced to avoid the divergence in the numerical calculations. As discussed in Ref. 0|, the dimension of k' is 
three and that of k, in the meson case is two. This extra dimension in k' should be caused by nonperturbative 
diagrams which include the form factor effects at low momentum region. We expect that k' ^ Aqc/jk, since 
^QCD is the only parameter related to confinement. In the numerical calculations, we let k' vary between 



0.01 GeV^ and 0.06 GeV^ [8 



39 



49 



The light quark propagator can be written as the following form 



lf]iM + qi - (jJi + le rjiM + qi + uji - lei 

where M and v are respectively the mass and velocity of the doubly heavy baryon, qi = q-v and q^ = q^ — qiv^ 
are the longitudinal and transverse projections of the relative momentum (q) along the doubly heavy baryon 
momentum (P), respectively, = y^m^"^^^, and A^^ are the projection operators given by 



. ---'- ^ (45) 

In the leading order of a 1/ mq expansion, the propagator of the scalar heavy diquark can be written as 

SD{q2) - , ' —-y (46) 

2mD(r]2M - qi - mo + le) 

After writing down the most general form for the BS wave function and taking into account its property 
under parity transformation, we can parametrize the BS wave function for the doubly heavy baryon with a 
scalar heavy diquark and a light quark in the following form: 

xp{q) = + 4tgs2)u{v), (47) 

where u{v) is the spinor of the doubly heavy baryon, and are the Lorentz-scalar functions of q^, qi and 
P2 = M\ 

Defining gs^^2^{qt) = /^5si(2)j O'^^ finds that the BS scalar wave functions satisfy the coupled integral 
equations as follows: 

~ , . f d^q't (m; + LOi)[Vi{qt ~ qj) + 2mDF{Q^)V2{qt - q't)] . , 
9si[qt) A 7nJ ^ 5siWtJ 

d'q't Vi{qt^q't)-2mnF{Q^)V2{qt-q't) , . , 

A 77? ^ 9t • qt 9s2[qt), (48) 

- , ^ / d'q't Vi{qt-q't) + 2mnF{Q^)V2{qt-q't) . , 

9s2{qt) = - / jTTT^ 7T7 ^ 9sAqt) 



d^q't imi~cui)[Viiqt-q't)-2mDF{Q^)V2iqt-q't)]qfq't - 



~9s2{q't)- (49) 



(27r)3 AuJimoiM - niD - c^i) qt 

The normalization condition for the doubly heavy baryon with a scalar heavy diquark and a light quark is 
given by (after imposing the covariant instantaneous approximation on the kernel) 



d^g d^g^ 



d 



Xp{q,s) Tr^Ip{q,q') xp{q' , s') ^ S^s' , (50) 

10 -To J 
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where ii(2) and ji(2) represent the color indices of the heavy diquark and the light quark, respectively, s^'^ is 
the spin index for the doubly heavy baryon and r^^^^^^'^ is the inverse of the four point propagator defined as 
follows: 

ir'''Hq,q') = S^'''S^''H2n)^S\q^q')Sr\qi)S^\q2). (51) 
After some algebra, Eq. (|50p can be written in the following form: 

-qJ -^^[T^'c[Xpiqt)Xp{qt)Siiqi){~'if]i^)Si{qi)SD{q2)] 

+Tr[xp{qt)Xp{lt)i~2tr,2)q2 ■ eSiiqi)SUq2)]} = 1, (52) 
where e = (1,0), Xrilt) and XpiQt) are the transverse projections of the BS wave functions given as follows: 

xp{qt) = -■iSi{qiy^xp{q)SD{q2y^, (53) 

and 

^p{qt) = ^tSD{q2)-'xp{q)Si{qi)-\ (54) 

respectively. 

Then, one can derive the transverse projections of the BS wave functions from Eqs. ([5^ (|54p . respectively: 



and 



where 



xpiqt) = [hsiiqt) + 4ths2{qt)]u{v), (55) 



Xp{qt) ^ u{v)[hs^{qt) + iths^iqt)], (56) 



and 



After substituting Eqs. (|55p ()56p into Eq. ()52p and integrating out the longitudinal momentum qi, the 
normalization condition can be written in the following form: 



hsMt) ^ I j^mqt-q',)+2mnF{Q')V2{qt-q^]gsAqt). (57) 



hsM) = I ^mqt-q',)-2mDF{Q')V2{qt-q't)]^~gsM)- (58) 



d\t 1 



:|(m; + uJi){2ri2Mujf + Ei]i{rnimD 



(27r)3 24Mm£,wf (-M + mo + cj;)^ 
-uJirriD + 2miUJi + Mui - Mmi) + Erii{-M + ttid + 2uji)ql)hl^{qt) 
+4ujfq^{r]2M + rnE)hs,{qt)hs,{qt) 
+ {2ri2M{mi ~ uji)qfiuf + Ei]i{-M + mo + 2uji)qf 

+Er]i{mi ~ u:i){~'mimD - ioimo - 2miUJi + M{mi + uji))q^hl^{qt)^ = 1, (59) 

where E = P ■ e. 
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Now let us define the BS wave function for the doubly heavy baryon composed of an axial- vector heavy 
diquark and a light quark as follows: 

x^(yi,2/2) = {omyMMlP) 

= (60) 

where ^^(2/2) stands for the axial- vector heavy diquark field. 

The BS equation for the doubly heavy baryon with an axial-vector heavy diquark and a light quark is given 

by 

X^(g) = Si{q,) I -0^G,.{P,q,q')xW)S'D'iQ2), (61) 

where S'^{q2) is the propagator of the axial- vector heavy diquark and Gp^ is the kernel for the BS equation 
given by 

- iGp^{P, q, q') = -QpJ ^ IVi{q, q') - 7" T^p,V2{q, q'), (62) 

where Tp,pv = [q^fx + (l2n)9pvP{Q^) is the effective vertex for the axial-vector heavy diquark coupling to the 
gluon, which was derived in Sec. III. 

In the leading order of a l/mg expansion, the propagator of the axial- vector heavy diquark can be written as 

'-^'^^ = -'2mn{mM-q,-mo + ^ey ^''^ 
Similar to the case of the doubly heavy baryon containing a scalar heavy diquark, we can parametrize the 
BS wave function for the doubly heavy baryon containing an axial- vector heavy quark and a light quark in the 
following form: 

Xp{q) = {gvi + 4t9v2)u''{v), (64) 

where u^{v) is the spinor of the heavy baryon, g„i and g„2 are the Lorentz-scalar functions of q^, qi and 
p2 = m"^. When the spin of the doubly heavy baryon is 1/2, u'^{v) = -^(t^ + v^^)'y^u{v), while the spin of the 
doubly heavy baryon is 3/2, u^^{v) is the Rarita-Schwinger vector spinor. 

After defining gvn2) (it) — I ^9vi(2) > can find that the BS scalar wave functions satisfy the coupled 
integral equations in the following: 

9M - - J iu;,moiM - run - CO,) '^'^"'^ 

~ f d'q't Vi{qt-q't) + 2mDF{Q^)V2{qt-q[) . , 
9v2\qt) — - TTTTT A 7U \ 9vi[qt) 

f d^ql {mi-uJi)[V,{qt-qi)-2mnF{Q^)V2{qt-ql)]qfq't -. .... 
J (27r)3 4uJimn{M-mn-u;i) qf ^''''> 

The normalization condition for the BS wave function for the doubly heavy baryon with an axial- vector heavy 
diquark and a light quark is given by 

^'^^ J (27r)8 Xp(g''^)[g^-^P^'-(g'g)J Xp{q,s) = 6,s', (67) 
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where Ip^^J^''^ is the inverse of the four point propagator defined as follows: 

ir;r'iq,q') = 5^^'^5^^'H27:)^6\q~ q')Sr\q^)S^l^{q2). (68) 
After some algebra, Eq. (I67p can be written in the following form: 

/ (^;;|l{Tr[XP^(<70lpp('7t)'5'i(<7i)(-«'?i^)'5i(gi)"5'2,''('Z2)] 
-Tr[xPa(gt)^Pp(9t)(«'72)(2(72 • eg,,!, - e^,q2y - g2M£i.)5z(gi)55''(g2)5^''(g2)] | = 1, (69) 
where XPa{qt) and XppiQt) are the transverse projections of the BS wave functions given by 

XPaiqt) - -tSi{q,r\Uq)SDc.a{q2)-\ (70) 

and 

Xppiqt) = -iSDp0{q2)~'^Xp{q)Si{qi)~'^ , (71) 

respectively. 

From Eqs. ([70)1 ([7T|) we get the expressions for the transverse projections of the BS wave functions as follows: 

XPa{qt) = [hvtiqt) + itK2{qt)]ua{v), (72) 

and 

^ppiqt) = up{v)\Ji^^{qt) + itK2{qt)], (73) 

where 

KAqt) = j ^mqt~q[)+2mDF{Q^)V2{qt-q[)]gM). (74) 

and 

= [ ^[Viiqt-q[)^2mnFiQ^)V2iqt~q[)]^~gM- (75) 

J (2t^) qt 

After substituting Eqs. (17^ ([751) into Eq. and integrating out the longitudinal momentum qi, the nor- 
malization condition can be written in the following form: 

/ (2^)3 24Mm^a;f (-M + run + {^"^' + -i)i^V2Mu:f + Er,,{m,mn 
-uJiTTiD + 2mi(jJi + Mui - Mmi) + Erii(-M + mo + 2uji)q^]hl^ {qt) 
+4ujfq^{r]2M + mE)KAqt)hvMt) 
+{2r]2M{mi - 0Ji)qtLjf + Et]i{-M + itid + 2uji)qf 

+Er]i{mi - uJi){-mimD - bJimo - 2m;W; + Af(m; + ijJi))qjhl^{qt)^ = 1. (76) 

In our calculation, we take the constituent masses of the light quarks as to„ = = 0.33 GeV , mg = 
0.45 GeV , and the scale of nonperturbative interaction Aqcd — 0.2 GeV. In order to solve the coupled integral 
equations p8|) (|49l) , we discretize the integration region into n pieces (with n sufficiently large) . In this way the 
integral equations are transformed into coupled matrix equations for the n dimensional vectors 5si(2) • Then, it 
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is easy to obtain the eigenvalue equations for ffs^j, ■ The same method is apphed in deahng with the coupled 
integral equations (|55|) After solving the eigenvalue equations, we obtain masses of the doubly heavy 

baryons shown in Table ITU The mass of the Sec obtained in our model is consistent with the experimental 
value, 3518.9 ± 0.9 MeV [50|. The obtained BS scalar wave functions for the doubly heavy baryons composed 
of a heavy diquark and a light quark are displayed in Fig. 4 and Fig. 5. One finds that the masses and the 
amplitudes of BS scalar wave functions for the doubly heavy baryons are independent of the spins of both the 
heavy diquarks and the doubly heavy baryons. In fact, this is just the consequence of the heavy diquark spin 
symmetry. So, we redefine gi = gs{v)i ^nd g2 = 9s(v)2 Fig- 4 and Fig. 5 for convenience. 



TABLE II: Values of the masses of baryons containing two heavy quarks. The lower (upper) masses correspond to 
k' = 0.01 (0.06) GeV. 





^bb 






^bb 


^bc 




M {GeV) 


10.08~10.10 


6.83~6.85 


3.52~3.56 


10.18~10.19 


6.94~6.95 


3.62~3.65 



V. Non-leptonic decay of doubly heavy baryons 

In this section, we will apply the obtained BS wave functions to calculate the non-leptonic decay widths for 
the doubly heavy baryons emitting a pseudo-scalar meson in the BS formalism. The Hamiltonian describing 
such decays reads [sif 

H,ff = ^V^bVJaiaiOi + 0202), (77) 

where Vcb and V^J^ are the elements of Cabibbo-Kobayashi-Maskawa matrix, U and D stand for the fields of the 
light quarks involved in the decay, Oi = [D-f'^{l ~ 7^)?/] [c7cr(l - 7^)6] and O2 = [c7°"(l - 7^)[/] [i57<j(l - 7^)6]. 
ai and 02 in Eq. (77) are defined as the linear combination of Wilson coeflScients (ci and C2), ai = ci + 02/ Nc 
and 02 = C2 -l- ci/Nc, where Nc is an effective number of colors which includes non-factorizable color-octet 
effects in the hadronization process. Due to the lack of knowledge about hadronization, ai and 02 are treated 



as free parameters and determined by fitting experimental data [52, l53|. Since & — >■ c decays are energetic, 
the factorization assumption can be applied in our calculation. Hence the decay amplitude of the two body 
non-leptonic decay becomes the product of two matrix elements: one is related to the decay constant of the 
factorized pseudo-scalar meson and the other is the weak transition matrix between the initial and final doubly 
heavy baryon states, 

(sg(P/)P(fc)|sg(PO) = ~^GFV,,Vi^a, Jd^ziPikmzhM-l'Wizme'^'' 

(sg(P/)|&(^)7'^(l-7')c(.)|sg(/^0), (78) 

where P{k) stands for the pseudo-scalar meson with momentum fc, Sgg, stands for the doubly heavy baryon 
composed of a (an) scalar (axial-vector) heavy diquark and a light u (or d) quark. 

The first matrix element on the right hand side of Eq. (|78l) is related to the decay constant of the pseudo-scalar 
meson, /p, which is defined as 



(P(fc)|5(z)7,(l - 7')t^(2)|0) - -ifpK. 



(79) 



15 



120 




|qj (GeV) 



FIG. 4: The normalized BS scalar wave functions for the doubly heavy baryons (Hqq/ ) composed of a heavy diquark 
and a light quark {u, d). The solid (dotted) lines are for k' = 0.01 (0.06) GeV. The upper solid and dotted lines are for 
Qiiqt) in unit of GeV~^. The lower solid and dotted lines are for gi{qt). 
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|qj (GeV) 



FIG. 5: The normalized BS scalar wave functions for the doubly heavy baryons (Qqqi) composed of a heavy diquark 
and a strange light quark (s). The solid (dotted) lines are for k' = 0.01 (0.06) GeV. The upper solid and dotted lines 
are for g2{qt) in unit of GeV~^. The lower solid and dotted lines are for gi(qt). 
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The decay amplitude (Eq. ([75)) ) is classified into three cases according to different initial and final doubly 
heavy baryon states as follows: 

Ml = {EQ,{Pf)P{k)\E*Q,{m, (80) 
M2 - {E*QAPf)P{k)\E*Q,{m, (81) 
Ms = {E*Q,{Pf)P{k)\EQ,{P,)). (82) 



p(k) 



c(p'i) 



(qi) 



FIG. 6: The Feynman diagram for the non-leptonic decays of doubly heavy baryons emitting a pseudo-scalar meson (P), 
taking the decay amplitude A4i for instance. EQc{Pf) and Sqi,{Pi) stand for the states of the final and initial doubly 
heavy baryons with momenta Pf and Pi, respectively. 4>Qc{Pdj) and ^QbiPoi) stand for the diquark states involved in 
the final and initial doubly heavy baryons with momenta Pd^ and Poi, respectively. &(pi), c(p'i), Q{p2), and l{qi) stand 
for different quark fields with corresponding momenta. P{k) stands for the emitted meson with momentum k. 



Let us first calculate the decay amplitude Mi. The Feynman diagram for Ali is shown in Fig. 6. A^i is 
related to the BS wave functions through the following equation: 

Ml = --^GpVcbVuDaifpk^ J d'^{xiX2yiy2ZiZ2z)e^'^'' 

XPj{x2-,xi)Si^{xi - yi)xp,(2/i,y2)5'^^(2;2 - ^^^^^^(za - ^2) 
{Q\Mzi)b{z)-i%l - 75)c(z)^^,(z2)|0), (83) 

where XPf ^-nd Xp stand for the BS wave functions for the final and initial doubly heavy baryons, respectively, 
(j)Qc and (pgj, are the diquark field operators involved in the final and initial doubly heavy baryons, respectively. 

The matrix element in Eq. ([83]) is related to the BS wave functions for the heavy diquarks through the 
following equation: 

{Q\Mzi)Kz)l%l - 7')c(^)^Qb(^2)|0) 



A 



Xg j ^^(271)4^4(^2 -p^)Tr[5(-p'J^^i^(p;)5(p2)^S^(p*)5(-pi)(l-7')7ir^'^ 
■^e-'-(P^i-Pn.+*^)e-^-^P^j e-^-2Pn^ ^ (84) 

where Xp^^ a-nd XPr!. a-re the BS wave functions for the heavy diquarks involved in the final and initial doubly 
heavy baryons, respectively, and mp). are the masses of the heavy diquarks involved in the final and initial 
doubly heavy baryons, respectively, Pof^^^i — Pd^^ ■ Vi {vi denotes the velocity of the initial doubly heavy 
baryon) and Pj^^^.^t = Pd^-) " PDj{i)i''^i the longitudinal and transverse projections of the final (initial) 
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heavy diquark momentum {Poffif ) along the initial doubly heavy baryon momentum (Pi), respectively, and ^'^(-^^ 
is the A-th polarization vector of the axial-vector heavy diquark involved in the initial doubly heavy baryon. 
Substituting Eq. (f84| into Eq. ([83]) . the decay amplitude A4i can be written in momentum space as follows: 

1 r A^nd'^n' 

Ml = -^GpV,,V^j,a^fpk,{2n)^d\P,- Pf ^ k) j ^-^{2^f5\q, - q[) 



X 



Analogously, we can derive the other two decay amplitudes (A^2 and A^s) in the BS formalism as follows: 
M2 = ~-j=GFV,,V{}j,arfpk„{2^f5^{P,-P}-k) j ^^{2nfS\q^ ~ q[) 

x{2nfSiPD,.i - ^f^^~Ph^t)5{PD.i - ^^^^^~Ph^t)XP,Sq't)Si{qi)xP.A<lt) 



4E / '^^{^^fs\p,-p',ms{-p',)i^^^^^^^^^ 



(86) 



and 



--J=GFV,,V^J,a^^pk,{27:f5\P, - ~ k) j -^(2nf5\q, - q[) 
x{2^f5{PD,i - ^ml^ - PI^MPd. - \lml^-PlMpfM't)Si{li)XpA'it) 



A' 



where ^'^'^^ ^ is the A'-th polarization vector of the axial- vector heavy diquark involved in the final doubly heavy 
baryon. 

The differential decay width for the two body decay reads 5C [ 



where D, denotes the solid angle, M stands for the decay amplitude. Mi is the mass of the initial baryon, and 
|k| is the absolute value of the three momentum of the particles in the final state in the rest frame of the initial 
state. 

Numerically, the parameters in the Hamiltonian {Gp, Vet, and T4/n), and the masses and the decay constants 



of the pseudo-scalar mesons are taken to have the following values (50|: Gp — 1.16637 x 10 ^ GeV ^, Vc 



Vcb 



0.0412, Vud = 0.97418, Vus = = 0.2255, Vcs = 0.9742, = 0.1396 GeV, tuk = 0.4937 GeV, mo = 
1.8696 GeV, mo, = 1.9685 GeV, U = 0.1304 GeV, fx = 0.1555 GeV, fo = 0.2058 GeV, and /d, 0.273 GeV. 
Our predictions for the non-lcptonic decay widths for the doubly heavy baryons emitting a pseudo-scalar meson 
are shown in Tables IIIIl IIV[ IVl and IVIl where the superscripts ^ and | denote the spin of the doubly heavy 
baryons. 
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TABLE III: Predictions for the non-leptonic decay widths for the doubly heavy baryons emitting vr meson. 

process decay width (10""a?Gel^) process decay width (10""a?GeV) 

r{E'J hI^tt) 0.343~0.362 r{n*J nl-jv) 0.591~0.607 

r{E*J ^E*Jn) 0.205~0.211 r{n*J ^ n*J n) 0.380~0.381 

r{Ell~^E'qn) 4.110-4.234 r(n;| ^ O;,^^) 7.606-7.643 

r(H* ^ E'Jtt) 0.848~1.101 T{n^^ Q,*jTr) 1.708-1.876 

r{E'J ^E*Jn) 0.415-0.587 riQ,*J ^n*jTv) 0.965-1.019 

r(S*| ^ eJItt) 8.626-12.110 T{n*J ^n*jTv) 19.435-20.529 



TABLE IV: Predictions for the non-leptonic decay widths for the doubly heavy baryons emitting K meson. 

process decay width (10"^^a?Ge\/) process decay width {IQ-^^alGeV) 

r{E*J EgK) 0.265-0.267 r{n*J -> ^t^) 0.469-0.482 

r{E'J~^E'jK) 0.165-0.171 r{n*J -^n'^K) 0.307-0.308 

riEll^Elfx) 3.317-3.425 r(n;| ^ fi^Jif) 6.150-6.170 

r(Hj^ -> E*JK) 0.649-0.845 r(fi4 ^ ^ll K) 1.313-1.557 

r{E*J -^E*Jk) 0.328-0.466 r{Q,*J ^ n'J K) 0.767-0.806 

r{E*J -^eIIk) 6.810-9.612 r{Q,*J ^ n'J K) 15.436-16.237 



VI. Summary and discussion 

In the heavy quark limit, a doubly heavy baryon can be regarded as a bound state composed of a heavy 
diquark and a light quark. We first establish the BS equations for both the heavy diquarks and the doubly heavy 
baryons, respectively, in the leading order of a l/mq expansion. The kernel for the BS equation contains the 
scalar confinement and one-gluon-exchange terms, which are motivated by the potential model and successfully 
used in the cases of mesons and heavy baryons containing a single heavy quark. Since the size of the heavy 
diquark is enhanced by In^mg with respect to l/mg, we also introduce a few form factors to the effective vertex 
for the heavy diquark coupling to the gluon in order to reflect the inner structure of the heavy diquark. 

The BS equations are solved numerically under the covariant instantaneous approximation, which is suitable 
for the weakly bound states of both the heavy diquark and the doubly heavy baryon. The obtained masses 



of the doubly heavy baryons are consistent with those from the lattice simulations Sj, 5^. It is found that 



TABLE V: Predictions for the non-leptonic decay widths for the doubly heavy baryons emitting D meson. 

process decay width {W-^^alGeV) process decay width (IQ-^^alGeV) 

r{E'J ^ ElD) 0.818-0.832 r{Q*J QlD) 1.404-1.466 

r{E'J^E'jD) 0.693-0.741 r{n*J^n*jD) 1.292-1.319 

r(S*(f ^ H*ji3) 13.885-14.834 r{Q,*J ^n*jD) 25.873-26.419 

r(H|^^H*|l>) 1.136-1.528 r{Ql^ ^ Q,*J D) 2.434-2.552 



r{El^ ^E'JD) 0.945-1.464 r{Q.l;? ^Ql^D) 2.383-2.426 

r{E*J -^E*Jd) 19.525-30.028 r{Q,*J ^n'jD) 47.895-48.741 
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TABLE VI: Predictions for the non-leptonic decay widths for the doubly heavy baryons emitting Ds meson. 



process 


decay width (10""a?GeV) 


process 


decay width (lO'^a^GeV) 






0.285~0.290 






0.511~0.515 


. 1 

r(sj - 




0.253~0.271 


mil^ - 


^^l^Ds) 


0.471~0.483 






5.065'-5.435 


ml! - 


-^^i!ds) 


9.437~9.668 


mi - 


> e'Jds) 


0.368^0.497 






0.794~0.828 






0.328~0.513 


ml! - 




0.828~0.851 


mil - 


• 3 

D ) 

' 1 — tec ^ s / 


6.763~10.520 


ml! ' 




16.650~17.099 



the properties of both the heavy diquarks and the doubly heavy baryons are independent of their spin in the 
leading order of a 1/ mg expansion. 

As we know, the superflavor symmetry relates doubly heavy baryons to heavy mesons, and hence the form 
factors of the transitions of doubly heavy baryons are reduced to the Isgur-Wise function, which is well known 



for heavy mesons [101, nil, 1561. The calculation of the doubly heavy baryon transitions is greatly simplified under 
the superflavor symmetry at the cost of ignoring the derivation from non-pointlike spatial dispersion of the 
heavy diquark. In this work, we directly calculate the decay amplitudes for the doubly heavy baryons using the 
BS wave functions obtained for both the heavy diquark and the doubly heavy baryons, instead of employing 
the superflavor symmetry. We give the predictions for the non-leptonic decay widths of doubly heavy baryons 
emitting a pseudo-scalar meson. Our results will be tested in the future experiments. 

In our calculation, for the propagators of the heavy diquarks and the light quarks involved in the bound states, 
we simply assume the forms of free propagators with the masses of the heavy diquarks and the light quarks 
taken to be the constituent ones. Actually, the real propagators should be solved using the Dyson-Schwinger 
equation. In such an approach, one has to guarantee the consistency between the kernel of the B S eq uation and 
that of the Dyson-Schwinger equation as required by the axial- vector Ward-Takahashi identity 2J]. This is a 
very complicated procedure and needs more careful investigations in the future. 

At the HERA-B and Tevatron facilities more than 10^ events involving double charm baryons are expected, 
while at the LHC one can expect about 10^ events 57|. Since the available energy at the LHC is much higher 
than the masses of Sec and Sbf,, it is believed that their production rates should be comparable. The decay 
widths of the doubly heavy baryons we present in this work will be tested in the forthcoming experiments. 

Since the BS equations are established at the leading order in a l/mg expansion, wc do not distinguish the 
different spins of both the heavy diquarks and the doubly heavy baryons. Such differences should happen at 
©(I/toq). l/mg corrections will be studied in the future. 



Acknowledgments 



One of the authors (M.-H. Weng) is very grateful to Dr. X.-H. Wu for a number of and very helpful discussions. 
This project is supported by the National Natural Science Foundation of China (Project Nos. 10675022 and 



21 



10975018) and the Special Grants from Beijing Normal University and by the Australian Research Council. 



[1] M. Mattson et al. (SELEX Collaboration), Phys. Rev. Lett. 89, 112001 (2002). 

[2] A. Ocherashvili et al. (SELEX Collaboration), Phys. Lett. B 628, 18 (2005). 

[3] R. Chistov et al. (Belle Collaboration), Phys. Rev. Lett. 97, 162001 (2006). 

[4] B. Aubert et al. (BABAR Collaboration), Phys. Rev. D 74, 011103 (2006). 

[5] M. Anselmino, P. KroU, and B. Pire, Z. Phys. C 36, 89 (1987) and references therein. 

[6] D. Ebert et al, Z. Phys. C 71, 329 (1996). 

[7] X.-H. Guo and T. Muta, Phys. Rev. D 54, 4629 (1996). 

[8] X.-H. Guo, A.W. Thomas, and A.G. Williams, Phys. Rev. D 59, 116007 (1999). 

[9] X.-H. Guo, K.-W. Wei, and X.-H. Wu, Phys. Rev. D 77, 036003 (2008). 

[10] M.J. White and M.,L Savage, Phys. Lett. B 271, 410 (1991) and references therein. 

[11] M.A. Sanchis-Lozano, Nucl. Phys. B271, 440 (1995). 

[12] X.-H. Guo, H.-Y. Jin, and X.-Q. Li, Phys. Rev. D 58, 114007 (1998). 

[13] S.-P. Tong et al, Phys. Rev. D 62, 054024 (2000). 

[14] W.-S. Dai, X.-H. Guo, H.-Y. Jin, and X.-Q. Li, Phys. Rev. D 62, 114026 (2000). 

[15] D. Ebert, R.N. Faustov, V.O. Galkin, and A. P. Martynenko, Phys. Rev. D 70, 014018 (2004). 

[16] J.M. Flynn and J. Nieves, Phys. Rev. D 76, 017502 (2007). 

[17] E. Hernandez, J. Nieves, and J.M. Verde- Velasco, Phys. Lett. B 663, 234 (2008). 

[18] W. Lucha, F. Schoberl, and D. Gromes, Phys. Rep. 200, 127 (1992). 

[19] E.E. Salpeter and H.A. Bcthe, Phys. Rev. 84, 1232 (1951). 

[20] N. Nakanishi, Suppl. Prog. Theor. Phys. 43, 1 (1969). 

[21] G.C. Wick, Phys. Rev. 96, 1124 (1954). 

[22] T. Nieuwenhuis and J.A. Tjon, Few Body Syst. 21, 167 (1996). 

[23] G.V. Efimov, Few Body Syst. 33, 199 (2003). 

[24] CD. Roberts, M.S. Bliagwat. A. H6U, and S.V. Wright, Eur. Phys. J. Special Topics 140, 53 (2007). 

[25] N. Nakanishi, Phys. Rev. 130, 1230 (1963). 

[26] N. Nakanishi, Graph Theory and Feynrnan Integrals (Gordon and Breach, NewYork, 1971). 

[27] K. Kusaka and A.G. WiUiams, Phys. Rev. D 51, 7026 (1995). 

[28] K. Kusaka, K. Simpson, and A.G. Williams, Phys. Rev. D 56, 5071 (1997). 

[29] V. Sauli and J. Adam Jr, Phys. Rev. D 67, 085007 (2003). 

[30] V. Sauh, J. Phys. G 35, 035005 (2008). 

[31] V.A. Karmanov and J. Carbonell, Eur. Phys. J. A 27, 1 (2006). 

[32] J. CarboneU and V.A. Karmanov, Eur. Phys. J. A 27, 11 (2006). 

[33] V.A. Karmanov and J. Carbonell, Nucl. Phys. B (Proc. Suppl.) 161, 123 (2006). 

[34] J. CarboneU, V.A. Karmanov, and M. Mangin-Brinet, Eur. Phys. J. A 39, 53 (2009). 

[35] J. Carbonell and V.A. Karmanov, Eur. Phys. J. A 46, 387 (2010). 

[36] J. CarboneU and V.A. Karmanov, hep-ph/1012.0246vl. 

[37] P. Jain and H.J. Munczek, Phys. Rev. D 48, 5403 (1993). 

[38] V.A. Miransky, LA. Shovkovy, and L.C.R. Wijewardhana, Phys. Rev. D 63, 056005 (2001). 

[39] X.-H. Guo and X.-H. Wu, Phys. Rev. D 76, 056004 (2007). 

[40] Z.-X. Xie, G.-Q. Feng, and X.-H. Guo, Phys. Rev. D 81, 036014 (2010). 

[41] H.-Y. Jin, C.-S. Huang, and Y.-B. Dai, Z. Phys. C 56, 707 (1992). 

[42] Y.-B. Dai, C.-S. Huang, and H.-Y. Jin, Z. Phys. C 60, 527 (1993). 



22 

[43] Y.-B. Dai, C.-S. Huang, and H.-Y. Jin, Phys. Lett. B 331, 174 (1994). 

[44] C.-H. Chang, J.-K. Chen, and G.-L. Wang, Commun. Theor. Phys. 44, 646 (2005). 

[45] E. Eichten, K. Gottfried, T. Kinoshita, K.D. Lane, and T.-M. Yan , Phys. Rev. D 17, 3090 (1978). 

[46] J.M. Cornwall, Phys. Rev. D 54, 6527 (1996). 

[47] D. Ebert et al, Z. Phys. C 76, 111 (1997). 

[48] S.S. Gershtein, V.V. Kiselev, A.K. Likhoded, and A.I. Onishchenko, Phys. Rev. D 62, 054021 (2000). 

[49] X.-H. Guo and H.-K. Wu, Phys. Lett. B 654, 97 (2007). 

[50] C. Amsler et al. (Particle Data Group), Phys. Lett. B 667, 1 (2008). 

[51] G. Buchalla, A.J. Buras, and M.E. Lautcnbacher, Rev. Mod. Phys. 68, 1125 (1996). 

[52] X.-H. Guo and A.W. Thomas, Phys. Rev. D 58, 096013 (1998). 

[53] O. Leitner, X.-H.Guo, and A.W. Thomas, Eur. Phys. J. C 31, 215 (2003). 

[54] R. Lewis, N. Mathur, and R.M. Woloshyn, Phys. Rev. D 64, 094509 (2001). 

[55] N. Mathur, R. Lewis, and R.M. Woloshyn, Phys. Rev. D 66, 014502 (2002). 

[56] H. Georgi and M.B. Wise, Phys. Lett. B 243, 279 (1990). 

[57] A.V. Berezhnoy, V.V. Kiselev, A.K. Likhoded, and A.I. Onishchenko, Phys. Rev. D 57, 4385 (1998). 



